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Abstract 

Junction temperature of a prototype SiC-based bipolar junction transistor (BJT) was estimated by 
using the base-emitter voltage (V B e) characteristic for thermometry. The V B e was measured as a function 
of the base current (I B ) at selected temperatures (T), all at a fixed collector current (I c ) and under very low 
duty cycle pulse conditions. Under such conditions, the average temperature of the chip was taken to be 
the same as that of the temperature-controlled case. At increased duty cycle such as to substantially heat 
the chip, but same I c pulse height, the chip temperature was identified by matching the V BE to the 
thermometry curves. From the measured average power, the chip-to-case thermal resistance could be 
estimated, giving a reasonable value. A tentative explanation for an observed bunching with increasing 
temperature of the calibration curves may relate to an increasing dopant atom ionization. A first-cut 
analysis, however, does not support this. 

Source of the Problem 

Intelligent use of a transistor in power management applications requires control of its junction 
temperature to some specified limits. The commonly specified thermal conductivity, junction to case, is 
not always available for packaged experimental or prototype devices obtained from R&D laboratories. 
Packaging obviously prevents the direct observation of the semiconductor chip. The sole recourse then is 
to find a terminal electrical property that can be exploited for thermometry. This paper reports an attempt 
to estimate the thermal resistance of a packaged SiC based BJT, from observations of electrical 
characteristics under low and high duty cycle power switching. 

Procedure Explained 

The regulation or switching of collector current, I c , to a load is the most direct way of setting the 
temperature rise of the chip above its case. And the case is assumed mounted to a temperature-controlled 
heat sink. The base-to-emitter junction voltage, V BE , depends on the base current, I B , and the junction 
temperature T. And so V BE (I B , T) and Ic are two elements on which the temperature rise thermometry can 
be based. But of course there are complications. 

A key realization is that external measurements at the base and emitter leads include bulk voltage 
drops in the chip added to the junction voltage. Ohmic drops due to the flow of Ic can be both significant 
as compared to the junction V BE as well as temperature dependent in doped SiC. In SiC, usually the 
dopant atoms are not fully ionized at room temperature, leading to a significant temperature dependence 
of bulk electrical resistivity as the temperature is raised. And there may be other temperature 
dependencies. Therefore our V BE thermometer must be calibrated taking I c into account. And this Ic must 
have the same pulse height as will be finally used in a higher duty cycle mode to heat the chip internally. 
At room and selected higher case temperatures, as set by an external heater, the external V BE is plotted 
against I B , in a very low duty cycle pulse mode. The Ic is set to a fixed pulse height that is representative 
of practical application of the particular BJT, by say a resistive load and a dc power supply. And the I B 
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pulses should be kept sufficiently high to ensure saturated switching of the BJT and just long enough to 
get steady V B e and I B readings. If properly done, internal heating of the chip then should be negligible. 

The so generated thermometry curves can be used to obtain an estimate of the junction-to-case 
thermal resistance. The duty cycle now has to be substantially increased to get a measurable increase of 
the steady state chip temperature above that of the heat-sunk case, while maintaining the same I c pulse 
height. C omparison of the V BE pulse height with that of the calibration curves then gives the temperature 
The thermal resistance then follows from the known case temperature and the average power loss in the 
chip shown by the displayed waveforms. 


Basic Experiment 

All thermometry work was done with variations of the circuit shown in figure 1. This circuit was built 
primarily for establishing a turn-off transient that is as rapid as the tum-on transient, but modifications 
worked well also for thermometry. 



Figure 1 . — This switching speed test circuit was used for 
thermometry calibration work. In the higher duty cycle, chip 
heating case, the negative pulser and its 50 Q load resistor 
were simply replaced by a dc power supply. Also at higher 
duty cycle, other resistors were varied to reduce their 
dissipation. The positive pulser provides the constant pulse 
height I B . 
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Temperature calibration data are shown in figure 2 below for a 4 A, 600 V rated SiC BJT made by 
United Silicon Carbide, Inc., and labeled P2. The decrease of V B e with increasing T at constant I B is 
equivalent to the usual shifting to the left of the I vs. V curves for say silicon P-N junction diodes, with 
the I-axis drawn vertically. But the curve bunching seen in figure 2 is not usual. However, in SiC devices 
the dopants are usually not fully ionized at room temperature. Hence as the temperature is raised, the 
degree of ionization increases and the body resistance drops. And the terminally measured V BE includes 
the resistive drop due to Ic, which thus decreases with increasing temperature until the dopants are fully 
ionized. The degree of ionization also plays a more subtle effect through the complicated physics of the 
so-called reverse saturation current I 0 in the diode equation. In fact, the simple diode equation, when 
solved for the junction voltage in the form 

V = (k T/q) ln[(I/I 0 ) + 1], (1) 

can be misleading as to temperature effects, because the complicated T-dependence of Io is not manifest. 
Being an aside, details of the mechanism that can (and usually do) generate a negative SV/ST in the diode 
equation are ferreted out in appendix A. The bunching of the curves seen in figure 2 is not, however, 
predicted there. This bunching is in any case undesirable, as it reduces temperature resolution and here 
especially in the region of interest around 200 °C. 

The P2 device was next mounted to a temperature controlled heat sink plate and set up to switch a 
pulsed 3 A current into a 100 Q non-inductive load resistor that was also mounted to the same heat si nk . 
Pulses of about 640 ns were spaced 5 ps, making a duty cycle of 0. 13. The combined heat from the load 



Figure 2. — Temperature dependence of the V BE versus I B curves. The duty cycle here is 
sufficiently low such that the chip and case temperatures are practically identical. 
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resistor and the transistor caused the temperature of the heat sink to slowly rise. When the Vbe pulse 
height matched the V B e of the 100 °C calibration curve (see fig. 2), waveforms were recorded, plate 
temperature (80 °C) was noted and operation ceased. Hence the case-to-chip temperature rise was then 
20 °C, 

Various scope measurements could then be done on the recorded traces. The exponential rise-time for 
the Ic was T on = 68.3 ns and the linearly approximated fall-time was x 0 ff = 134 ns. The formulas giving the 
energy loss in these approximations are 

E on = (1/2) V 0 Io Ton (exponential approximation) 


and 


E 0 ff = (1/6) Vo Io T 0 ff (linear approximation), 

where V 0 and I 0 (= I c ) are the switched voltage and current. The arithmetic then gives E on = 3.16x] 0 7 J 
and E 0 ff = 2.07xl(T 5 J. According to the scope measurement, the on-state Vce was only 0.4 V, whereas 
according to the static characteristics of the P2, about 4 V was expected at 100 °C and 3 A. Using the 
scope measurement, the conduction loss is only 7.6x1 CT 7 J, which is quite small compared to the 
switching loss. The base-drive energy was also small, being about 0.15x10 J. Finally, adding all losses 
and dividing by the 5 ps period gives an average power of 1 1 W into the chip. 

The junction-to-case thermal resistance R eJ c is thus about 20/1 1 = 1.8 °C/W. This is somewhat lower 
than that of a roughly comparable silicon BJT, such as the planar BUH50 (50 W, 800 V, 4 A). But such a 
comparison can only be very rough, because the P2 has no given power rating. Nevertheless, the value of 
the thermal resistance obtained appears to be reasonable and gives some confidence in the method. 

Final Remarks 

The initial experimental evidence supports the feasibility of the presented method of temperature rise 
thermometry. But it also points out the problem of curve bunching and consequent loss of temperature 
resolution seen in figure 2. The root cause in the described experiment appears to be the increase in body 
conductivity as more dopants in SiC (especially the acceptors) become ionized with increasing 
temperature. This problem may be a characteristic of a material, such as SiC, whose commonly used 
dopant atoms are not fully ionized at room temperature. However, a workaround may be possible by 
using pulse control methods.. 


Alternative Schemes 

In a modified procedure, consider generating the low duty cycle, calibration curve set (fig. 2.) with 
the collector lead open. Making I c = 0 hopefully should reduce the nonlinear bunching with increasing T 
of the V B e versus I B curves, if dopant ionization is at the heart of the matter. The general idea is to use 
higher duty cycle pulses of I c to substantially heat the chip, but measure a steady V BB and I B while I c = 0. 
For example, after a sufficient number of base drive pulses to cause a substantial temperature rise of the 
chip, the collector power supply could be effectively disconnected, for the duration of one base pulse or 
one period. This would be a matter of pulse control logic. A scheme might utilize a commercially 
available high voltage pulser for the collector supply, with base drive applied continuously or pulsed in a 
synchronism. Values of Vbe and I B could be measured during the Ic dead time. In all such schemes, a fast- 
recovery diode (GaAs, Schottky) is needed in the collector lead to block current flow in the base-to- 
collector direction. 


NASA/CR— 2006-2 14258 


4 



Conduction and Switching Losses 

Occasionally the concern has been voiced that higher on-state conduction losses in a SiC based, as 
compared to a Si based, BJTs would be seriously adverse to their application as switches. The present 
example of the P2 device, however, shows that this is circumstantial and not always so. Therefore, let us 
look a little closer at the comparison of conduction losses to switching losses expected for a waveform of 
essentially square pulses of width w, repeated with a period x. The duty cycle D then is D = w/x. Let the 
switched voltage and current be V 0 and Io. 

Approximating both tum-on and turn-off by the same exponential function, the total switching loss in 
one period then is 


Es — Vq Io x s , 


where the exponential switching time constant x s is related to the usual scope measured 10 to 90 percent 
switching time x s , 10 - 9 o by x s = x s ,io- 9 o/hi 9. The conduction loss in one period is simply 

E c = I 0 “ RcE.on W , 

where RcE.on is the collector-to emitter on-resistance, mimicking MOSFET terminology. Division of E s 
and E c by x then gives the average switching and conduction loss powers Ps and Pc- The ratio of these 
powers thus is 


Pc/P s -E c /E s -I 0 Rce,o„w/(Vo x s ) . 

To get a numerical feel, assume a switching frequency of 200 kHz and a duty cycle of 0.5, giving 
w = 2.5 ps. The x s ,i 0 - 9 o is not likely to be significantly less than 50 ns even for SiC power BJTs, giving x 
s = 23 ns. At least for the P2 device, the RcE.on was observed to be about 1.5 D. at room temperature on a 
static curve tracer, but appeared to be less under fast pulse conditions. 3 A is about the rated current for 
this device, but operation at 300 V is at about ‘A of the rated voltage. With these values, Pc/Ps = 1-6. But 
operation at higher frequency (smaller w) and higher voltage is clearly possible for this device and 
V, io- 9 o = 50 ns is likely to be on the short side. Hence the evidence shown here does not support the claim 
of necessarily unacceptably high conduction losses in SiC based BJTs. 
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Appendix A 

Temperature Dependences in the Diode Equation Applied to SiC 

The usual observed decrease of V with increasing T at constant I means that in the simplest p-n 
junction diode modeling equation, i.e., Eq. (1), the so-called reverse saturation current I 0 must have a 
T-dependence of at least exponential order. We shall review here this dependence, pointing out the source 
of the usually negative (5V/5T)|i and the effect of dopant degree of ionization. A caveat is that Eq. (1) 
may be inadequate if there is significant recombination in the depletion region, in which case additional 
terms are required. 

A standard form for I 0 , written in terms of the electron and hole diffusion constants, minority carrier 
diffusion lengths, intrinsic carrier density and the ionized dopant densities, is 


Iq — q A 


Dn n? ^ Dp n? 
l n n a LpNJ, 


(Al) 


The above Io can be put in more fundamental terms such as the temperature, electron and hole 
mobilities, densities and effective masses, and minority carrier lifetimes and the band gap energy. Use is 
made of the Einstein relation and standard textbook formulas [1] for n ; , etc., to finally give 

Io =4Aq 1/2 h -6 (2 tt) 3 (m* m*) 3/2 


En 


\ 1/2 


+ n _ 


Ep 

X D 

V P J 


T/2 


(kT ) 7/2 e kT 


(A2) 


where A is the junction area, q is the magnitude of the electronic charge, p = N^ and n = N jt, . In the case 

of an NPN transistor, the emitter region will be more heavily doped than the base and so the n in the 
bracket can be neglected. 

Some of the terms in the above expressions can vary widely, even in Si. The lifetimes are sensitive to 
temperature and the presence of recombination centers. The mobilities vary as T”“, where 2<a<2.5. The 
effective masses are not particularly temperature sensitive. Unlike in Si, the dopants used in SiC are far 
from being fully ionized at room temperature and hence their degree of ionization must be calculated 
from their positions in the bandgap and the temperature. The following values have been suggested [2]: 

m n * = 0.5 m 0 x n = 100 ns p n = 900 cn^wV 1 N A « 3xl0 17 cm 3 

m p * = 0.3 m 0 T p =100ns p p = 90 cn^VV 1 J « 500 A/cm 2 

E c - E d = 0.05 eV E A - E v = 0.25 eV 


In spite of uncertainties, a rough estimate for SiC of the terms in Eq. (Al) or Eq. (A2) can be obtained. 
Let us write Io in the form 


E, 


I 0 =Af(T)e"KT , 


(A3) 


where the function f is defined by 
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f(T)=4q l/2 h 6 (27r) 3 (nin nip) 3/ “ 


(A4) 


/ \ 1/2 

Tn_ 

l X „J 


-n 


-l 


f V/2 

Mp_ 

v t p. 


(kT) 


7/2 


and the mobilities, lifetimes and carrier densities are functions of T. 

The most temperature-sensitive terms in Eq. (A4) are assumed to be the p and n, as the degree of 
ionization involves an exponential dependence on T. Finding only the p is sufficient here, since the 
emitter is assumed to be heavily doped and the donors sit much shallower in energy. The ionization 
formula is developed in reference 3, the result for acceptors being 



/ 

1 + 

V 


M/2 


4Na 
N £ 


-1 , 


(A5) 


where 


N ^=-^-( 2 7rm; KT/ h 2 ) 3/2 e“A^ , (A6) 

and the degeneracy factor g D is usually set to 2. If N A »N^, such when nearing carrier freeze-out, then 
p « (N a N c ) 1/2 is a good approximation. Thus at 300 K, where = 1.39xl0 14 , only p = 6.46x1 0 15 cm 3 
out of the 3xl0 17 cm -3 acceptors are ionized, compared to p = 5.80xl0 16 cnf 3 at 500 K. The 
corresponding change in body resistivity is obviously quite large. 

For any reasonable forward current, the “1” in Eq. (1) can be neglected and hence Eq. (1) can be 
written in the form 


qV = E g + k T ln[J/f(T)] , (A7) 

where J = I/A is the current density and E g is the band gap (about 3 eV in SiC). From this follows the 
sensitivity 


(0V/5T)|j = (K/q) [ln(J/f) - (T/f) (df/dT)] , (A8) 

of the junction voltage to temperature. 

A rough estimate in SiC shows that at room temperature, the value of J/f(T) is considerably below 
unity, making ln[J/f(T)] negative. Indeed, f(300) « 1.4xl0 7 A/cm 2 , giving (K/q)ln(J/f) ~ -8.8x1 (T 4 V/K, 
which is much below the roughly -1 .3x1 0 2 V/K observed. Moreover, the uncertainties in mobilities, 
lifetimes and ionized dopant densities seem insufficient to correct for this discrepancy. 

Let us ferret out more systematically the main T-dependences hidden in the function f and 
contributing via the df/dT. The simplest such in Eq. (A4) is some power of T. Besides the manifest T 7/2 , 
the mobilities too have been modeled as proportional to T~ 2 ' 3 . The contribution to (<3V/0T)|j in Eq. ( A8) 
from total power law T“ proportionality in f is 

-ff ( a(fa) a( g T a) )=- -f a=-8.61xl(Ha (V/K) . (A9) 
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But from Eq. (A4) and the above remarks, a ~ 1 and surely below 10. Hence this contribution is rather 
small and constant. 

The exponential T-dependence of p seems likely to contribute more. Its complications are best ironed 
out step by step. After some tedious algebra, one finds that 


dp _ i SN^ 

dT~ 2 ST 


Nr 


. 1/2 


-1 


1 + 4 


N a 


Nr 


(A10) 


and also that 


SN^ _ i 

ST _ T 


( 3 , AE a a 
v 2 kT , 


N r 


(All) 


where AE a s E a - E v . Again in the case of N A »N A which is applicable here at 300 K, 


Sp 3 AE a 
ST ~ 2T I 2 + kT 


(n A nJ /2 4 


^3 , AE a ^ 

v 2 kT j 


(A 12) 


Using the above results, the contribution to (SV/ST)|j arising from Sp/ST is found to be 


kt 

Sf Sp^ 


( \ 

kT 

'f Sp_" 


( 3 , AE A 

qf 

v SpST y 


qf 

VC / 

[pSTj 

~ 2q 

U kT j 


(A13) 


within the approximation based on N A »N^. At 300 K, the value of this positive contribution is 4.79x10 
V/K. 

For general interest, the above contributions can be totaled. Omitting considerable algebra, the final 
form is 


-4 


/ ^ \ 


r r 

8V 

K 

^-a + ln 

W 

u" q l 

4 L 


Co(kT) 


,11/4 


(A 14) 


where 


Co- 2(27r)9/4 gDq l/2 h -9/2 Nx l/2( I] 


*\3/2/ A3/4 


A. V /2 


(A 15) 


For the given numerical values (and g D = 2), C 0 = 1.28xl0 65 A m 1 (Joule) 11/4 . At K = 300, and 
J = 500 A/cm 2 , (SV/ST)|j = 8.61 xlO 5 [3/4 - a +ln(4.36xl0~ 3 )] V/K = 8.61 xlO 5 (3/4 - a -5.44) V/K. 

This diode equation based analysis applied to SiC predicts that (SV/ST)|j is, with high confidence, 
negative and increasing in magnitude with increasing temperature at room temperature. It provides no 
support for the curve bunching seen in figure 2. 
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Worth a passing note is the absence of AE a in Eq. (A14). This is actually apparent from the outset 
from Eqs. (A7), (A4) and (A6), whenever p ~ (N A NJ 12 is a valid approximation. The icTln operation in 

ae a 

Eq.(A7) then reduces the e 2kT factor, arising from Eq.(A6), to -AE a / 2, which contributes zero to 

E g 

(oV/5T)|j. The band gap factor e kT appearing in Eq.(A2) suffers a similar fate. The fact that E g and AE a 
may themselves have a temperature dependence was simply ignored here. 
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